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Abstract 

Color-octet contributions to quarkonia production at LEP are studied herein. The 
short distance formation of heavy quark-antiquark pairs in color-octet configurations via 
gluon fragmentation processes is significantly enhanced relative to the creation of color- 
singlet pairs via heavy quark fragmentation. But the subsequent long distance hadroniza- 
tion of these colored pairs into physical quarkonium bound states is suppressed compared 
to the nonperturbative evolution of their colorless counterparts. We find that the overall 
LEP rates for gluon fragmentation into prompt Upsilon and Psi vector bosons exceed those 
from heavy quark fragmentation. Inclusion of the dominant color-octet quarkonium pro- 
duction channel eliminates sizable discrepancies between previous predictions and recent 
measurements of prompt Z ^ J /ij} -\- X , Z ^ + X and Z —> T + X branching fractions. 
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1. Introduction 



Rare decays of Z bosons to heavy quarkonia final states are currently under study at 
LEP . The predominant majority of charmonia observed at the electron-positron facil- 
ity originate from B meson decays. A smaller fraction come from prompt processes which 
create charm-anticharm pairs at the primary interaction point. These two qualitatively 
different sources of cc mesons can be distinguished via vertex detection measurements. 
In the 66 sector, no delayed production mechanism exists, for bottomonia are among the 
heaviest hadrons created in Z decays. Bottomonia consequently result only from prompt 
reactions. 

During the past few years, a number of theoretical predictions for quarkonia pro- 
duction rates at LEP have appeared in the literature. The heavy quark fragmentation 
processes 

Z^Q^Q Z^Q+Q 

I and I _ (1.1) 



have been claimed to represent the main source of prompt Upsilon and Psi vector bosons . 
The rates for these modes are orders of magnitude larger than those for non-fragmentation 
channels such as the gluon radiation reaction Z i^qgg- But the recent first detection 
of Upsilons by the OPAL collaboration casts doubt upon the assertion that the decay 
processes in ( |1.1|) account for most bottomonia produced at LEP [jl|. In particular, the 
measured branching ratio 

3 

^ Br(Z ^ T{nS) + X) = (l.2+[j-| ± 0.2) x 10"^ (1.2) 

n=l 

is ten times larger than the prediction based upon bottom quark fragmentation. Similarly, 
the numbers of prompt J/t/j and t/j' mesons which OPAL has observed |0] 

Br(Z ^ J/i; + X)„, = ^t^f'"^^"' x Br(Z ^ J/^ + X)totai 

3 n -i_ o Q^ XX 1 n— 4 



(0.077 ± 0.058) X (3.9 ± 0.2 ± 0.3) x 10"^ = (3.0 ± 2.3) x 10 



A^(V'')prompt 



Br(Z ^^' + X)prompt = r^'^ X Br(Z ^ V' + ^)totai 

VP J total 



(0.14 ± 0.09) X (1.6 ± 0.3 ± 0.2) x 10"^ = (2.2 ± 1.5) x 10"^ 

(1.3) 
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exceed the abundance expected from the modes in (|1 . Ij ) . Although these OPAL findings 
are prehminary and have yet to be confirmed by other LEP groups, the large discrepancies 
between the data and existing theoretical estimates are interesting. 

In this article, we investigate a potential resolution to this puzzle. We suggest that the 
largest contribution to prompt quarkonia production stems from gluon rather than heavy 
quark fragmentation. This possibility was first considered in 1991 by Hagiwara, Martin 
and Stirling 0. The branching ratio calculated by these authors for the gluonic process 

Z q + q+g* 

I (1-4) 

turned out to be smaller than the branching fraction found later for the quark mode in 
( p. . 1|) . But their result was based upon the assumption that the entire reaction in ( |1.4|) takes 
place at an energy scale greater than or equal to the i/'q mass. Since Z — > qqg* — > qqtpqgg 
decay starts at 0(q;^), it is not surprising that the rate for the purely short distance part 
of this transition is small. 

Theoretical understanding of gluon fragmentation into heavy quarkonia has signifi- 
cantly matured since the time Hagiwara et al. carried out their analysis. High energy for- 
mation of heavy quark-antiquark pairs in color-octet configurations followed by low energy 
evolution into color-singlet hadrons is now known to represent the dominant source of large 
p±_ quarkonia at hadron colliders IfzHT^ . Inclusion of gluon fragmentation mechanisms into 
charmonia and bottomonia cross section calculations appears to eliminate orders of mag- 
nitude discrepancies between theoretical predictions and recent Tevatron measurements 
1^,14|. However, it is important to test this new understanding of quarkonia production 



in other settings. We are thus motivated to study the implications of color-octet mechanism 
ideas for Z decays to Upsilon and Psi final states. 

Our paper is organized as follows. We first calculate the LEP rates for charmonia and 
bottomonia production arising from lowest order color-octet processes in section 2. We 
then review the heavy quark fragmentation prediction for both color-singlet and color-octet 
quarkonia production in section 3. We examine next in section 4 the gluon fragmentation 
contribution io Z ^ T X and Z ij) ^ X. Finally, we compare all these different 
production modes in section 5 and close with a summary of our findings. 



2. Lowest order color-octet quarkonia production 

Quarkonia production at LEP involves several different length scales. The first is set 
by the Z boson mass Mz which represents the e^e~ collider's current operating energy. 
The second is fixed by the mass Mq of the quark or antiquark inside a quarkonium state. 
The heavy constituents are bound together over a distance scale determined by their 
momentum MqV, and they interact with each other over time periods set by their kinetic 
energy MqV^ . Since the velocity v ~ 1/ log Mq of the Q and Q inside the bound state 
is small compared to the speed of light, these scales are all widely separated from one 
another: 

{MQv'^f < [M^vf <^Ml<^ Ml (2.1) 

The presence of so many different energy scales complicates the analysis of quarkonium 
production. But the most important characteristics may be simply distilled within an effec- 
tive field theory framework called Nonrelativistic Quantum Chromodynamics (NRQCD) 



which systematically keeps track of the hierarchy in (2T) [15|. This effective theory is 
based upon a double power series expansion in the short distance strong interaction fine 
structure constant and the small velocity parameter v. We will work throughout this 
paper within the NRQCD framework. 

To begin, we consider the decay mode 

Z ^g + QQ[^s+^Lf] (2.2) 

which creates a colored heavy quark-antiquark pair recoiling against a hard gluon. The 
lowest order diagrams that mediate this transition are illustrated in fig. |1|. The spin, 
orbital and total angular momentum quantum numbers of the QQ pair are indicated in 
spectroscopic notation inside the square brackets in ( |2.2|) , and its color is labeled by the 
octet superscript. Using the projection techniques described in refs. |lT3| , |T6| , pr7[| , we can 
start with the general Z ^ g + Q + Q amplitude and pick out terms which correspond to 
the formation of a QQ pair in a specified partial wave, spin and color state. For example, 
the short distance production amplitude for a colored L = 0, 5" = 1 pair is given by 



A{Z{k + P) ^ ga{k)+QQ[^sf\{P)) 



short 

2QnQr, M (2-3) 
3 



Here M = 2Mq represents the pair's mass, §2 and gs denote the SU{2),^ and SU{'i)c 
gauge couphngs, T3 stands for the third weak isospin eigenvalue and 6w equals the weak 
mixing angle. The outgoing Q and Q in fig. |] propagate nearly on-shell with a combined 
momentum P. But the intermediate heavy quark or antiquark line which attaches to the 
decaying Z is off-shell by an amount of 0(M|). As we shall see, the intermediate particle's 
extreme virtuality strongly suppresses the lowest order color-octet process in ( |2.2| ) relative 
to quark and gluon fragmentation reactions that take place at higher order in the dg 
expansion. 

After the colored QQ'\^^^'^L^j^] pair is created within a spacetime interval set by 
the Z boson mass, it evolves over a much longer time and distance scale into a physical 
quarkonium bound state. The pair transforms into a colorless hadron via the emission 
or absorption of one or more soft gluons. These long wavelength gluons bleed off the 
color of the QQf^^+^Lj^-*] object, but they carry away only 0{Mqv'^) worth of energy and 
momentum [|15[. In the Mq 00 limit, this hadronization process should be perturbatively 
computable. However in the real world where Mqv'^ ~ Agcij, the amplitude for a heavy 
color-octet pair to evolve into a physical quarkonium state is nonperturbative. The long 
distance matrix element therefore cannot be computed from first principles and must 
instead be determined from experiment. 

Combining together the square of the short distance amplitude in ( |2.3| ) with the long 
distance factor MsO^q) = \-^{QQ[^S[^^] 'i/'Q(n)) ^^^^p, we obtain the partial width 

nz 9 + - +X) = ^ f 1 - ^) (2.4) 

I y^vvL 1 J ^^Q^ ) 3 gin2 cos2 0^ V 2 )\ Ml J M, ^ ' 

which is averaged (summed) over initial (final) polarizations and colors. The color-octet 
squared amplitude decomposes into an infinite sum over NRQCD matrix elements 

+ E {0\O^'''^'^\'S,)\0)Br{xUm)^Mnh) + ---} 

J=0,l,2 

(2.5) 

which determine the probabilities for a QQf'^S'j^^''] pair to turn into a i/jq through all possible 
intermediate channels. Numerical values for these NRQCD matrix elements have recently 
been extracted from CDF Upsilon and Psi cross section measurements in refs. [lT3| , p^ . 
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Once the universal matrix elements are known from one experiment, they may be applied 
to the study of color-octet quarkonia production at any other. 
Inserting the long distance squared amplitude values 

Ms{J/^) = 6.8 X 10"^ GeV^ 
MM') = 2.0 X 10"^ GeV^ 

(2.6) 

3 ^ ' 

^ A^8(T(n^)) = 6.4 X 10"^ GeV^ 

n=l 

into (|]D along with the parameters = 1.48 GeV, Mb = 4.88 GeV, as(2Mj = 0.28 
and asi2Mb) = 0.19, we find the Z ^ J/tp + X , Z ^ ^j' + X , and Z ^ J2n '^(nS) + X 
branching fractions listed in the first row of Table I. Comparing these predictions with the 
OPAL numbers in eqns. ( p..2|) and (|1.3|) , we see that the contribution to Upsilon and Psi 
production at LEP from the color-octet process in ( p.2|) with an L = 0, = 1 intermediate 
state is negligible. Yields from similar reactions involving colored QQ pairs with different 
orbital and spin angular momenta should be even smaller since their hadronization into 
physical bound states takes place at higher orders in the velocity expansion. We therefore 
conclude that 0{as) color-octet quarkonia production channels can be ignored without 
loss. 



3. Heavy quark fragmentation 

Parton fragmentation represents the dominant source of heavy quarkonia in high en- 



ergy collisions [|T8[ . As we have seen, the propagators of short distance intermediate partons 
are usually off-shell by 0(l/i?^) where E denotes some characteristic energy of the hard 
reaction. But in fragmentation processes, the quark or gluon which splits into a quarko- 



nium state plus other partons is typically off-shell by only 0{1/M^). At sufficiently high 
energies, fragmentation contributions to quarkonia cross sections dominate over all non- 
fragmentation competitors since the former are enhanced relative to the latter by powers 
of E'^/Mg. This outcome is guaranteed to occur even if the short distance part of the frag- 
mentation reaction takes place at higher order in the as expansion, for additional factors 
of the QCD coupling only logarithmically suppress reaction rates with increasing E. 

One of the first applications of these fragmentation ideas was to charmonia and bot- 
tomonia production at LEP. In ref. 0, Braaten, Cheung and Yuan calculated the prob- 
abilities for charm and bottom quarks to respectively split into Psi and Upsilon bound 
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states. These authors then folded together the fragmentation probabihties with the Z ^ cc 
and Z ^ bb partial widths and determined the prompt quarkonia branching fraction 
Br(Z — > i/jq + ^)Qfrag- We will review the findings of Braaten et al. below and then 
generalize their color-singlet results to include color-octet contributions as well. 

Heavy quark fragmentation production of Upsilon and Psi bosons at LEP proceeds 
at 0{a1) through the Feynman graphs illustrated in fig. |^. S The mass s of the virtual 
heavy quark or antiquark which splits into the outgoing QQ pair ranges over the interval 
3Mq < s < Mz — Mq. ipQ production is maximized when the pair is created near the Mq 
scale in a color-singlet L = 0, = 1 configuration. In the fragmentation approximation, 
the contribution of this channel to the Z ^ ipQ partial width simply equals the product of 
the Z QQ decay rate and the hadronization probabilities Dq^^^ and Dq^^^ ||]: 

r(Z ^ Q + Q + QQfs['^] + X) = T{Z ^ QQ) X {Dq^^^ + D^^^^) (3.1) 



where 



30 



57 In 2 



M2 



(3.2) 



All nonperturbative information associated with -i/^g bound state formation resides within 
the squared amplitude 



5: {(0|O^^-)(^^i)|0)Br(^,(m)^^,(n)) +■..}. ^^"^^ 



m>n 



Numerical values for the color-singlet NRQCD matrix elements which enter at lowest order 
in the velocity expansion into this long distance factor can be extracted from the charmonia 
and bottomonia wavefunctions at the origin tabulated in ref. fl^ . Adopting the values 



Mi{J/tp) = 6.1 X 10-2 GeV^ 

Miit/;') = 2.9 X 10-2 GeV^ 

^ ^ 3.4 

3 

^7Wi(T(n>S)) = 2.3 x 10-^ GeV^, 



^ The two diagrams displayed in fig. |2| dominate in axial gauge over two other Z — > 'ipgQQ 
graphs which do not possess an obvious fragmentation interpretation. Analysis of ipq production 
via heavy quark fragmentation is most transparent in this particular gauge. 
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we derive the heavy quark fragmentation branching ratios hsted in the second row of 
Table 1. 

It is instructive to compare these color-singlet results with their color-octet analogues. 
The partial width for Z Q + Q + QQ[^s[^^] i/jq + X is readily obtained from 
r{Z — > Q + Q + QQ[^s[^^] — s> i/^Q + X) by performing the following color- factor and 
squared amplitude substitutions in eqn. (|3.2|) : 

(^^) and Mii^^) ^ MsO^q). (3.5) 

These alterations diminish both the short and long distance contributions to the color-octet 
rate relative to its color-singlet counterpart. The color-octet heavy quark fragmentation 
branching ratios appearing in the third row of Table I are consequently much smaller than 
their color-singlet counterparts. 

Comparing the predicted rates for the processes in eqn. (|1.1|) with the Upsilon and 
Psi measurements in eqns. ( |1 . 2| ) and ( p.. 3D , we see that heavy quark fragmentation does 
not appear to account for most prompt quarkonia observed at LEP. Although the data 
are preliminary and their error bars are large, the sizable mismatch between theory and 
experiment suggests that other rapid production mechanisms are at work. We therefore 
turn to consider gluon fragmentation contributions to Upsilon and Psi production in the 
following section. 



4. Gluon fragmentation 

A number of striking disagreements between theoretical predictions and experimen- 
tal measurements of quarkonia cross sections at hadron colliders have arisen during the 
past few years. The CDF collaboration has discovered that prompt J/i(j and i/;' produc- 
tion at the Tevatron exceed color-singlet differential rate predictions by orders of magni- 



tude [pO| , |2T[| . CDF has also observed greater numbers of T(15'), T{2S) and T(35') vector 
mesons than were originally anticipated [^. These experimental surprises from Fermilab 
have stimulated several recent theoretical investigations of quarkonia creation in general 
and gluon fragmentation processes in particular [0-0. The incorporation of color-octet 
production mechanisms into cross section calculations has been demonstrated to yield sig- 
nificantly better descriptions of charmonia and bottomonia measurements than those based 
upon color-singlet mechanisms alone. Since numerical values for long distance color-octet 
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matrix elements are a priori unknown, the absolute magnitude for da{pp ^ i^Q + X)/dp± 
must simply be fitted to data. But the shapes of the predicted and measured differential 
distributions at high p± agree quite well once gluon fragmentation reactions are taken into 
account. Moreover, the fitted values for the color-octet matrix elements are consistent with 
NRQCD velocity counting rules. Color-octet contributions thus appear to eliminate large 
discrepancies between theory and experiment. 

Given that gluon fragmentation plays a central role in quarkonia production at the 
Tevatron, it is interesting to study its impact at LEP as well. Gluon fragmentation me- 
diated decays of Z bosons into Upsilon and Psi final states start at O(a^) through the 
reaction displayed in eqn. (|1.4| ) and illustrated in fig. ^. A straightforward computation 
yields the differential rate for this process: 

+ 5^ [8EiE2{Ei + E2) - 2{3Ef + AE1E2 + 3E|) + 4(^1 + E2) - l] 
+ 5\l - 2i?i)(l - 2i?2)}/[(l - 2i?i)2(l - 2E2f] . 

(4.1) 

Here Ei and E2 represent the dimensionless lab frame energies of the quark and antiquark 
scaled relative to the Z mass, S = M/Mz equals the rescaled quarkonium mass and K 
denotes a constant prefactor 



= (4.2) 

which contains the rescaled long distance matrix element A^8('*/'q)/-^|- The 0(5^) and 
0((5^) terms inside the curly brackets in ( [4.1| ) are negligibly small compared to the leading 
0{6^) term since 5^ <C 1. Similarly, corrections to the Z —>■ qqg* partial width from the 
masses of the q and q are unimportant and have been ignored. 

Integrating the differential expression in ( [4.1| ) over the quark and antiquark energy 
limits < < (1 - d^)/2 and (1 - d^)/2 - Ei < E2 < [1 - d^/{l - 2Ei)]/2, we obtain 
the partial decay rate 

T{Z^q + q + g* ^^^+X) = ^{{l + S^f\ogSHog j^^^ + (3 + 46^ + ^5^) log 

+ 5(1 -n + 2(1 + 5^)2^- / 5' \ / 1 



(4.3) 
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The 1 jS^ term appearing alongside the constant prefactor K originates from the propagator 
of the virtual gluon. We again recall that the fragmenting g* is off-shell by an amount of 
order Mg rather than M|. The integrated rate for the gluonic process is thus enhanced 
by the same 1/5^ factor as its heavy quark fragmentation analogue. 

The gluonic channel's rate is further amplified by a large double logarithm. The source 
of the (log (5^)^ factor appearing in the first term of (|4.3|) can be traced to the double pole 
residing in the first term of (|4.1| ). The double log in the integrated rate results from the 
overlap of an infrared singularity with a coUinear divergence. As log 5^ is numerically 
large in both the charmonia and bottomonia sectors, higher order corrections to the gluon 
fragmentation partial width are likely to be significant. The leading logarithms should be 
resummed in order to obtain a more reliable estimate for the rate at which -i/^g quarkonia 
are generated via gluon fragmentation at LEP. Performing such a resummation is rather 
subtle, and we will return to this issue in a future work. For now, we will simply evaluate 
the lowest order expression with the understanding that subleading corrections may not 
be negligible. 

With the extra double log enhancement factor, the short distance contribution to the 
gluon fragmentation process in ( p..4|) dominates over that for the heavy quark fragmenta- 
tion reaction in ( |1 . 1|) . But the long distance color-octet matrix element A^sIV'q) fo^^ the 
former is numerically two orders of magnitude smaller than the corresponding color-singlet 
squared amplitude M.\{^^q) for the latter. The relative importance of these two different 
fragmentation mechanisms thus depends upon the outcome of a competition between short 
and long distance factors. To determine the winner, we evaluate the partial width in ( |4.3| ) 
and list the resulting Z ^ T and Z ^ i/; branching fractions in the last row of Table I. 
Comparing these entries with the others in the table, we see that the gluon fragmentation 
rate beats those for all other prompt quarkonia production channels. 

Inclusion of the dominant color-octet contributions brings theoretical branching frac- 



tion predictions into line with the experimental values quoted in eqns. (|1.2| ) and ( |1.3| ). 
These results consequently represent a nontrivial success of the color-octet picture. We 
also note that the rate for gluon fragmentation production of J/V' is consistent with the L3 
upper bound Br{Z qqg* J/t/j + X) < 7.0 x 10"^ [H and the more recent DELPHI 
limit Br(Z qqg* ^ J /ij) -\- X) < 4.1 x 10"'* 0. Improved LEP measurements of this 
branching ratio along with its ip' and T analogues would provide useful cross checks on 
NRQCD color-octet matrix element values. Angular distribution and spin alignment data 
would further test color-octet mechanism ideas. In short, further observation of Z decays 
to Upsilon and Psi final states would enhance understanding of basic quarkonium physics. 
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5. Conclusions 

In this article, we have investigated color-octet contributions to prompt quarkonia 
production at LEP. i The quantitative results of our study compiled in Table I can be 
qualitatively understood by considering the order of magnitude estimates displayed in 
Table 11. In the first column, we list the perturbative QCD order for each of the production 
channels that we have considered in this paper. We next recall in the second column the 
short distance factors associated with each mode. As we have seen, the 1/5^ = M|/M^ 
enhancement of the fragmentation reactions overwhelms their 0{as) suppression relative to 
lowest order non-fragmentation processes. Rates for the former consequently swamp those 
for the latter. The branching fraction for the gluonic mode in ( p..4|) is further amplified by 
a large double log. In the third and fourth columns of Table II, we enumerate the relevant 
color and flavor factors for each of the production channels. It is important to note that 
these factors enhance the gluon fragmentation mode by an order of magnitude relative to 
its heavy quark counterpart. Finally, we list the long distance matrix elements associated 
with each process. An interesting interplay between the various factors appearing in the 
columns of this table determines the relative importance of the production mechanisms 
listed in the separate rows. 

The nonperturbative suppression of color-octet processes generally implies that they 
are subdominant compared to color-singlet reactions. Intermediate 66['^5'^^^] and cc['^S'P'*] 
contributions to bottom and charm quark fragmentation into Upsilon and Psi final states 
are negligible compared to those from 66['^S'(^'*] and cc[^S'(^'*] pairs. But in certain circum- 
stances, short distance enhancement of color-octet modes can offset their long distance 
suppression. A striking example of this phenomenon has recently been proposed as a res- 
olution to the CDF -0' surplus problem [|12[. Short distance enhancement of g —>■ cc[^s[^'^] 
could yield 50 times more high p± t/j' charmonia than g cc['^S[^^] even though the 
hadronization rate for colored pairs from the first reaction is suppressed by 0{v^) com- 
pared to that for colorless pairs from the second. A similar phenomenon may play an 
important role at LEP. We have found that inclusion of the dominant gluon fragmen- 
tation channel (|1.4| ) removes sizable disparities between previous theoretical predictions 
and recent experimental measurements of Z ^ J/i(j + X, Z + X and Z ^ T + X 

branching fractions. Color-octet mechanism ideas thus appear to solve puzzles that have 
been encountered in both hadron and lepton collider settings. Our confidence is therefore 



^ Similar work has recently been reported by Cheung, Keung and Yuan in ref. |24]. 
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bolstered that these same ideas can be apphed to quarkonia problems in other contexts as 
well. 
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r TO H n p i" 1 on 

Mechanism 


Yir(Z — > 7/1/;") 




3 




1.6 X 10-^ 


4.6 X 10-* 


1.0 X 10-6 


^ _^ n -L n -L 00\'^si^^] 

ZJ 7 ^ ^ Ko^Ko^^ i^-^ J 


7 7 X 10~^ 


3 7 X If)-^ 


1.6 X 10-^ 


Z^Q + g + QQp5f^] 


1.1 X 10-^ 


3.2 X 10-8 


5.7 X 10-8 




3.3 X 10-^ 


9.6 X 10-^ 


4.1 X 10-^ 



Table I 



Production 
Mechanism 


QCD 
Order 


Short Distance 
Factor 


Color 
Factor 


Flavor 
Factor 


Long Distance 
Matrix Element 


Z 




as 


1 


2 


1 




z 


Q + Q + QQ[^S[^^] 


al 


1 

52 


V J 


1 




z 


Q + Q + QQ[^S[^^] 




1 


iV^-l 


1 




z - 


^q + q + QWs['^] 


Oil 


(log 52)2 
^2 


2 


5 





Table II 
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Figure Captions 

Fig. 1. Feynman diagrams which mediate the lowest order color-octet process 

Fig. 2. Feynman diagrams which mediate the heavy quark fragmentation process 
Z ^ Q + Q + QQ^^+'^L^P] ^ Vq + ^- Two other non-fragmentation graphs 
whose 0{a1) contributions are subleading in axial gauge are not pictured. 

Fig. 3. Feynman diagrams which mediate the gluon fragmentation process Z Q+Q+g* 
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